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Abstract 

We consider a 3-brane embedded in a warped 5-dimensional background with a 
dilaton and a Kalb-Ramond 2-form. We show that it is possible to find static solutions 
of the form of charged dS/AdS-like black hole with horizon which could have a negative 
mass parameter. The motion of the 3-branc in this bulk generates an effective 4- 
dimensional bouncing cosmology induced by the negative dark radiation term. This 
model avoids the instability that arises for bouncing brane in a Reissner-Nordstr0m- 
AdS bulk. 



1 Introduction 



Braneworld models [1,2] have generated, during the past decade, enormous attention, due 
to the dramatic change they inspired in our understanding of extra dimensions. According 
to this framework, our universe is a "brane" embedded in a higher-dimensional space, on 
which the Standard Model fields are confined, while gravity is localized near the brane 
by the warped geometry of the extra dimension. It is possible to construct models in 
which the brane evolution mimics a Friedmann-Robertson- Walker (FRW) cosmology, with 
modifications at small scales due to the gravitational effect of the bulk spacetime on the 
brane [3-5]. In particular, provided the bulk is taken to be a Reissner-Nordstr0m-AdS black 
hole, such modifications can lead to bouncing 4D cosmological models [6]. Unfortunately 
the brane, during its evolution in the bulk, always crosses the Cauchy horizon of the AdS 
black hole, which is unstable [7,8]. 

In this paper we present a different model, in which this problem is avoided. We consider 
a brane embedded in a supergravity background in which both the dilaton and the Kalb- 
Ramond 2-form are turned on (but without a dilaton potential). By dualizing the 2-form, 
we obtain Einstein-Maxwell like equations of motion, but with a different sign for the 
kinetic term of the Maxwell-like field. The static solution is therefore different, and the 
term that dominates at high curvature is like "stiff matter" with positive energy density. 
Even though this implies that the energy contribution at high curvatures is positive, so that 
it can not drive a bounce, this opens an interesting possibility of having negative energy 
contributions at intermediate curvatures, by letting the mass of the black hole be negative. 
The parameter space allows this while avoiding a naked singularity. In this case, we show 
that it is possible that the brane bounces before crossing the black hole horizon, so that the 
effective 4-dimensional cosmological evolution will not suffer the instability of [8]. 

The paper is organized as follows: In section [2] we derive the bulk equations and present 
the static solution, which is an asymptotically (A)dS black hole, charged under the Kalb- 
Ramond field with the dilaton frozen at its VEV. Section O is devoted to the study of the 
position of the horizons, and it is shown that in a certain region of the parameter space 
(in particular for a negative cosmological constant) there is no naked singularity even if 
the mass parameter of the black hole is negative. Then, in section H] we embed the brane 
and allow it to move, to mimic a cosmological evolution for a brane observer [9]. The 
occurrence of the bounce and its position is discussed here. Finally, in section [5] we present 
our conclusions and discuss open problems. 
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2 The 5-D solution 



Consider the low-energy string effective action: 
M3 



-Ts / dHV^e'^ (1) 



(with Habc = c^[yi-S_BC])i which describes a 3-brane embedded in a 5-dimensional bulk with 
dilaton and Kalb-Ramond 2-form [10,11]. We set the dilaton potential to zero for simplicity, 
and allow the dilaton to be non-minimally coupled with the cosmological constant, the 
antisymmetric tensor and the brane [12,13]. However, the brane is taken to be neutral with 
respect to the antisymmetric field. The presence of the brane action in ([1]) gives a singular 
part, which we will take into account by the Israel junction condition [9,14], so now we will 
consider only the bulk part. The equations of motion are: 



Gab 



2 ^ ^ 24 



gAB + dAcl^dB^ + ^c'^^'^'HaceHb^^, (2) 



VaV^0 - aie'^'^L - ^e'^^^H^ = 0, (3) 



Vc (e^'^^H^^^) = 0, (4) 

where H"^ = HabcH^^^. 

The equation for the antisymmetric tensor ^ can be solved by the ansat 

Substituting this into eq. (H]) we get 

Vc[e^^^,CABDE^^^^^-.,^^^^^ABCDE^^^^^^ = 1^^^^^^ (V^V^ - Vz.Vc) ^i. 

= y^^^^Rcni'AM = ^. (6) 

The last equality follows from the second Bianchi identity R^abc] ^ = 0- Now the equation 
for Am can be obtained by the identity f?[^-?^BCD] = 0) while we have to substitute the 
expression we get for H in eqs. ([2]) and ^ to get the correct equations of motion. Observing 
that 

H'^ = 12e~2'^2^2^ 
HamnHb^''' = 2e-2-2(2FVi?-4FAci^s^), (7) 



^See also [15], though our solutions differ from theirs. 
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where ^ 

Fmn = 2 C^mAn — VnAm) , (8) 



we get the foUowing equations: 



Gab 



9AB + dA<t>dB(t> - 2e-''^'t'FAcFB ^, (9) 



VaV^(^ - aie^^'f'L - -^e-^^^F"^ = 0, (10) 



Vb (e'^^'^F'^^^ = 0, (11) 



which look like the equations of motion of an Einstein-Maxwell model, but in this case the 
signs of the Maxwell fields are reversed. 

Now, following [9] we search for a static solution with a maximally symmetric 3-space. 
The metric ansatz is: 

ds^ = -m)dt' + + r^dn^) + (12) 

We also assume that the gauge field is purely electric, i.e. Am = (A(i?),O,0), and that the 
dilaton also depends only on the radial coordinate, (f) = <p{R). Then the Einstein equations 
read: 

2R R^ R^ 4 2-'^' 

f" f f k 1 1 

2 R R^ R^ 4 2^^ ^ ' 

It is easy to sec that, summing up the first and second equations, we get that the dilaton 
must be constant, <p'{R) = 0, so that we can reabsorb it into the definition of the coupling 
constant. So, we are left with only the vector field; the equations are 

^" + 34 = 0. (14) 
R 

As usual, by the conservation of the energy-momentum tensor, only two of these equations 
are independent. A solution of these equation can be cast in the form 

/(i?) = -^-JL + k--R\ 
■'^ ' 3i?4 i?2 6 

A{R) = ±-|. (15) 
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Notice that this solution is quite similar to the one that was found in [6], but the term 
proportional to the charge is opposite in sign. Because of this, it is possible to avoid a 
naked singularity even if fi is negative. In the next section we will discuss the location of 
horizons for different values of the physical constants. 



3 Horizons in the Kalb-Ramond Black Hole 

We can track the location of the horizons by finding the zeros of the metric function / in 
(jlSp. For simplicity we will assume the spatial part of the metric to be flat, k = The 
zeroes of the function / can be found by solving the equation: 



x^ + 6^x + 2- 



0, 



(16) 



with X = B? (so that we are interested only in positive solutions). This equation is conve- 
niently solved by use of the Chebyshev radicals. We can identify three different cases (let 
us stress that we will assume < from now on): 



Case 1 When 



L > 



we have one real solution 



X 



- cosh 
ij 



- cosh ^ 
3 



(-2^)3/2 J 



(17) 



(18) 



but this solution is always negative, so it is not acceptable. Therefore in this part of 
the parameter space the background has a naked singularity. 



Case 2 When 



< L < 



8^ 



we have three real solutions: 
Xl = 2 



2fi 



■ cos 



X2 
X3 



_2,/-|^cos 



-Xl - X2, 



arccos 



— arccos 
3 



_QWh_\ 

(-2^)3/2 J 
(-2^)3/2 J 



(19) 



(20) 



of which the second one is negative, X2 < thus unacceptable, and xi > x^. In this 
case we have two horizons, just as in the usual Reissner-Nordstr0m black hole. The 
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horizons can be written as: 



R 




Case 3 This case occurs when 

L < 0, 

and we have, again, one real solution, which can be written as 



(21) 
(22) 

(23) 



X 



X 



V 

2W ^ cosh 
V i-' 



1 

- arccos 
3 



(-2/^)3/0 J 



for ^ < L < 0, 

for L<^. (24) 



This solution is always positive, so that in this case we have only one horizon, which 
can be written as 

1/4 I 



'C 



2gV5L\ 



where the Chebyshev polynomial C1/3 is intended to be 



Ci/3(i) ^ 



2 cos (i arccos(|)) < t < 2 
2cosh (|cosh-^(|)) t>2. 



(25) 



(26) 



Since the bulk is asymptotically AdS, this is the most interesting case. As we will 
see in the next section, it is in this case that the cosmological evolution on the brane 
undergoes a viable bounce. 



4 Cosmological evolution of the moving brane 

In this section we embed a 3-brane in the bulk described previously, by cutting out R > a, 
and imposing a Z2 symmetry at the edge. We then let the brane move through the bulk, 
a = a{T), where r is proper time on the brane. The movement of the brane in the 5D bulk 
induces a cosmological evolution on the brane via the Israel junction condition [9]. To see 
this, wc need to calculate the extrinsic curvature of the brane. The unit vectors tangent 
and normal to the moving brane are: 
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UA = \a,0,-^^ I, (27) 



/ 

so that the (spatial part of) the extrinsic curvature is 



1 V7+ 



2 



Kij = lij, (28) 

where 'jf^u is the induced metric on the brane, which is FRW with scale factor a, and the dot 
represents a derivative with respect to proper time on the brane. The modified Friedman 
equations are obtained from the junction conditions, 

1 1. 



K^J = [T^j - -T^^,J . (29) 

where T^j_u is the brane energy-momentum tensor. For T^y = X'^^u + {p + p)u^Uy + pj^u, 
where / is the brane tension, we find that: 

rr2_ P ,1^4 ^ , ( P \\ ^ ^ C^O) 

where 

M4 = -, L4 = + -. (31) 

^ 3 M6 ' ^ 3 M6 2 ^ ^ 

li w = p/p is constant, then by energy conservation on the brane, p = poa^^^^^'^\ and the 

effective Friedmann equation (j30p can be expressed in terms of the sole a. 

The eff'ective Friedman equation contains, apart from the quadratic term in the brane 

energy density, a dark radiation term proportional to p [3,4] and a term that behaves like 

"stiff matter", proportional to Q [6]. But, as pointed out in the previous section, in our 

case this last term, which dominates at high curvature, has a positive contribution, hence 

it cannot drive a bounce in the cosmological evolution. Nevertheless, we can allow for a 

negative contribution from the dark radiation term, by letting p be negative. There is 

always an "intermediate" curvature regime in which the negative dark radiation dominates, 

so that it can be responsible for a possible bounce. Notice that the bounce would happen at 

a larger radial coordinate in the bulk than in the case discussed in [6] , so that the instabilities 

which seem to rule out the bounce [8] could be avoided. In order to discuss these issues, 

let us simplify the setup by considering a pure tension, spatially fiat brane. The Friedman 

equation reduces to 

^-Y + ^ + 3^- ^^^^ 

Notice that, the 4D cosmological constant becomes dominant at late time, so that, in order 
to have an asymptotically de Sitter universe, it must be positive. This is always true in the 
dS bulk, while in the AdS case the tension has to satisfy the inequality 



l>\l-hM^. (33) 
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A bounce occurs if the scale factor a{t) reaches a minimum: d = 0, a > 0. The first 
condition can only be met when fj, < 0. The behaviour of as a function of a is depicted 
in Fig. [T] for different values of the parameter 

Q^VlI _ Rkr ^2^^ 



2(-/u)3/2 U ' 

which corresponds to the ratio between the 4D de Sitter curvature radius and the charac- 
teristic length of the Kalb-Ramond black hole obtained by the charge to mass ratio. 





Figure 1: The Hubble parameter H as a function of the scale factor a. The three colors 
represent decreasing values of 15.81 (red, dotted), 1.41 (blue, dashed), 0.70 (black, 

solid. 



There are two values of a for which H = 0, but, of course, the only one compatible with 
a late time de Sitter evolution is the one located at the largest value of a. As the scale 
factor shrinks, the Hubble parameter follow the negative branch of the plot backwards to 
zero, then it start to grow again, following the positive branch. Note that this heuristic 
observations prove that the point at which the Hubble parameter goes to zero (for the 
largest branch) is actually a minimum of the scale factoi§|. Quantitatively, the zeroes of the 

^ Quite interestingly, the smallest branch that starts from a = can mimic the evolution of a closed 
universe, in which the scale factor reaches a maximum and then starts to decrease, even if the space geometry 
of the brane is assumed to be flat. But we will not discuss this issue any further. 
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equation ([32]) can be obtained by solving the equation: 

f 



x-^ + -^x + ^ = 0. 



(35) 



We require L4 > to have a well behaved late time evolution and fi < for the bounce to 
occur. We find that there are three real solutions to eq. (j35p when 



< L4 < 



Q 



4 ^ Rkr < h, 



(36) 



The solutions can again be expressed in terms of the Chebyshev radicals given by eq. ()20p 
with L replaced by 2L4. The second solution is negative so it has to be discarded. We are 
only interested in the larger of the two solutions xi. So we conclude that a bounce can 
actually occur at the scale factor 



V2 



ii 



1/4 



' COS 



1 



■ arccos 



■2(-/x)3/2; 



(37) 



Next we need to discuss the sign of the second derivative of a, to prove that is actually 
a minimum as expected. Taking the derivative of (j32p we find: 



a 
a 



2Q^ H L4 

--— ~4 ~' 



3 a6 a4 ■ 3 

The sign of o at the bounce can be deduced by studying the function: 



(38) 



3a^. 



3/x 2 2Q2 



12 I -A 



3/2 



cos 



— arccos 
3 



6^ 2 
L4 



'2(-^)3/2; 



3Q^ 

3g2 



(39) 



where the second equality comes from using ()35p . Now, observing that cos[. . .] > 1/2 and 
using inequality ([36|) we get 



L4 



> 2 



uj 



3/2 



> 2 



(L4)3/2 U 



0, 



(40) 



which proves that there is actually a bounce at a = at- 

The crucial requirement for a bounce to be acceptable is that it occurs before the brane 
crosses the black hole horizon. So we have to check that there exist a non-empty region of 
the parameter space in which ab > R+ or ab > Rq- 
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Figure 2: Plot of Oil) clS cl function of the tension I for different values of : 0.5 (black, 
darker), 0.02 (red, lighter). Dashed curves of the same color represent the position of the 
horizons for the same value of The case with 2 horizons is on the left, the case with 

one horizon is on the right. 

We discuss briefly the case in which the bulk black hole has two horizons, i.e. case 2 
of the previous section, so that the physical parameters satisfy the inequality ([19]). As we 
will see, this case has to be discarded. In fact, Fig. [2] (left) shows a plot of ah as a function 
of the brane tension and the position of the two horizons i?+ and R-. We see that the 
position of the bounce is always inside the external horizon, and coincides with it for the 
tensionless brane. Analytically, it is easy to see that is a decreasing function of I, and, 
since L4 ^ L/2 as Z ^ 0, we readily get at R+ in the tensionless limit. So, this case is 
expected to have an instability similar to the one discussed in [8] (even though in that case 
the brane actually bounces after crossing both horizons). 

Next, we discuss the most interesting case of a negative bulk cosmological constant 
and a single horizon, case 3 of the previous section. Fig. [2] (right) shows ab{l), and the 
corresponding value of Rq. In contrast to what happens with two horizons, there is a 
region, when the brane tension is close to the critical value expressed in eq. (j33p . in which 
the bounce radius is greater than the horizon position, so that the entire evolution of the 
brane lies in the physically viable region outside the horizon. This feature is quite general, 
and the reason is easy to understand, since we can see from ([37j) that Of, ^ 00 as L4 ^ 0. 
Analytically, we have to solve the inequality af, > Rq- After some algebra, we find that L4 
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has to satisfy the following inequalities: 




L4 < --CjJ2 



2^)3/2 I (-2/i)3/2 



) 



for 



for 



(-2/i)3/2 



(-2/i)3/2 



<5, 



>5, 



(41) 



(42) 



where the Chebyshev polynomials have been defined in (j26p . So there is always an allowed 
value of the brane tension for which the brane evolution lies entirely outside the horizon. 

Finally, to get the full dynamical evolution from a 4D perspective, we need to solve 
the effective Raychaudhuri equation psp . It can not be solved analytically, but numerical 
solution confirms that the scale factor undergoes a bounce from a contracting phase to an 
expanding one, as illustrated in Fig. [31 



Figure 3: Plot of the numerical solution of a{t). The initial conditions used for the numerical 
integration are a{t = —5) = 125, a{t = —5) = —125. 



5 Conclusions and outlook 

In this paper we presented a braneworld model in which the cosmological evolution of the 
brane is non-singular, and the brane universe bounces smoothly from a phase of contraction 
to a subsequent expanding phase. The cosmological evolution on the brane is induced by its 
movement through a static bulk AdS black hole supported by a non-trivial Kalb-Ramond 
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antisymmetric 2-form. The solution is similar to the Reissner-Nordstr0m-AdS solution 
presented previously in the literature [6], but in our model the bounce is induced by a 
negative dark radiation term sourced by the black hole mass term. In order for this to 
be possible, the integration constant fi which appears in the metric solution (jlSp . and 
which is proportional to the mass of the black hole, must be negative. Black holes with 
negative mass have been considered in the literature [16,17]. This does not lead to a naked 
singularity, as we have shown in section [3J But it may create some difficulties related with 
the overall stability of the solution (and possibly to its physical interpretation). Stability 
analysis of charged black holes in 5 dimensions has proven to be a very difficult task, 
and there appears to be no definitive answer for the standard Reissner-Nordstr0m-(A)dS 
black hole [18]. More recently, it has been proved [19] that a class of topological black 
holes with negative mass is stable in every dimension, which can be regarded as hint of 
the good behaviour of 5-dimensonal negative mass black holes against perturbations. The 
thermodynamical interpretation of the black hole solution requires the temperature of the 
black hole to be positive. The temperature can be related to the value of the derivative of / 
at the horizon, so the latter must be positive in order to have an acceptable thermodynamical 
behaviour. In our solution this is always the case, whatever the value of /i. Stability and 
thermodynamics of the bulk solution presented in this paper are worth further investigation, 
which we plan to do in a forthcoming paper. 

Future investigation also should extend the study to more general settings. For example, 
we have only considered the case of a spatially flat, pure tension brane. While the curvature 
term is not expected to change dramatically the picture, it is possible that the presence of 
matter could spoil the bouncing behaviour. In fact, the bounce is driven by the negative 
dark radiation term, so a sufficiently large amount of positive "ordinary" radiation could 
possibly compensate for the negative energy term, thus making the singularity appear again. 
The study of the case of a non-empty brane is also important to test the model against 
the observational constraints on the early universe evolution, such as nucleosynthesis. In 
practice, it is a nice feature that the bounce occurs at a large scale, but it should not be 
"too large" , so that the universe would have enough time to undergo its standard evolution. 

Another interesting development would be to consider a charged brane. The string 
theory embedding of the braneworld scenario would require the interpretation of the brane 
as a D-brane on which the open strings, which represent gauge particles, are confined. Of 
course, a full string theory derivation of braneworld models is still lacking. Nevertheless, a 
D-brane is naturally coupled to the Kalb-Ramond 2-form via, for example, a Dirac-Born- 
Infeld action, which can be implemented also in the present model. But in this case the 
ansatz ([5]) is not expected to hold, so the calculation becomes much more complicated. 

Bouncing cosmologies have been proposed in various contexts [20-33] as an alternative 
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to inflation. They have of course the remarkable feature of avoiding the initial singularity, 
but also provide for an alternative, but still perfectly acceptable, solution of the kinematic 
problems of the standard cosmological model, without invoking an unknown inflaton. Nev- 
ertheless, scalar fluctuations observed by the WMAP and SDSS [34] experiments seem to 
favour a nearly scale-invariant spectrum, which is in agreement with the prediction from 
chaotic inflation, but very difficult to obtain in a bouncing model. So, as a second step, it 
would be essential to study the behaviour of the scalar (and tensor) perturbations. 
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